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' Abstract. We study the Landau-Zener transition with the quantum Zeno effect 



in an open dissipative system populated by a large number of bosons. Given the 
quantum Zeno effect is strong enough, both discrete and continuous quantum Zeno 
measurements are found to stabilize the Landau-Zener transition. Both the a^-type 
longitudinal relaxation and cr z -type transverse relaxation in the bosonic system are 
analyzed as a model of continuous quantum Zeno measurements. While both of them 
improve the signal-to-noise ratio in terms of the ground state population, the er^-type 
relaxation can further boost measurement sensitivity and thus lead to a polynomial 
speedup with the number of bosons in the system. For a system that contains a large 
number of bosons such as in a Bose-Einstein condensate with more than 10 bosons, 
this equates to several orders of magnitude speedup. 
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1. Introduction 

In the quantum Zeno effect, if a quantum system is frequently projected in a certain 
measurement basis including the initial state, the system is pinned on this initial statepQ. 
Moreover, if the sequence of measurements is performed along a specific path, the 
system will evolve along this arbitrarily chosen path. This effect is called the inverse 
quantum Zeno effect [2j [3] and has been proposed to manipulate quantum systems 
by Kitanojl]. The measurement basis of the sequence typically follows the diagonal 
basis of the temporal Hamiltonian of an adiabatic evolution. The goal of the control 
sequence is to use the inverse quantum Zeno effect to keep the system in the ground 
state which is one of the eigenstates of the temporal Hamiltonian. In this sense, the 
inverse quantum Zeno effect can also solve the problems targeted by adiabatic quantum 
computation [5J. However, the physics behind the measurement-based evolution is a 
dephasing and dissipation with respect to the reservoir, while adiabatic evolution is 
completely coherent. Thus the two schemes yield different time complexity [3]. 

More recent developments have revealed that quantum computation by engineered 
dissipation and decoherence is a promising alternative to coherent quantum computation 
[HI El |8]. It is also well-known that bosonic systems are able to enhance the dissipation 
process. One example is the property of final state stimulation, whereby in a dissipative 
process, if N bosons are present in the final state, then the cooling rate is enhanced by 
a factor of N + 1. A classical computational scheme that takes advantage of bosonic 
final state stimulation was proposed in Ref. [9]. It was found there that the time for 
solving to a typical Ising model type problem could be accelerated by a factor of 1/N, 
where N is the number of bosons per Ising model site. 

In this paper we study the performance of quantum computation by the quantum 
Zeno effect in an open dissipative system for a many-boson system. Note that henceforth 
we do not discriminate between the quantum Zeno effect and the inverse quantum 
Zeno effect. We focus on the scheme involving measurements along the path of the 
Hamiltonian for a Landau-Zener transition |10| ITT], the simplest example of an adiabatic 
quantum computation. We show that the Landau-Zener transition can be stabilized 
by a series of strong enough quantum Zeno measurements. More importantly, the 
quantum Zeno effect is enhanced significantly in the presence of many bosons and thus 
can accelerate the scheme of the Landau-Zener transition by improving the measurement 
sensitivity. 

2. Discrete model 

First, we show that the Landau-Zener transition can be realized by the discrete 
quantum Zeno effect by a series of projective measurements. In the bosonic system, 
the Hamiltonian for the Landau-Zener transition is written as 

H = -e{t)j z - AJ X , (1) 
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where -T/2 < t < T/2, and e(t) = 2e t/T = vt. J z = {a) a - tib)/2 and 
J x = (a)b + b Jf a)/2, where a and b are the bosonic annihilation operators for the two 
energy levels, for example, up-spin and down-spin states, respectively. Suppose that 
e 3> A > 0, and the system of N bosons is initialized in the initial ground state 
|0) a |iV)fe. By evolving the Hamiltonian slower than the adiabatic limit, the system is 
pinned in the ground state and reaches the state \N) a \0) b at t = T/2. 

The same goal can be achieved by a discrete quantum Zeno effect, consisting of 
a series of projective measurements along the diagonal basis of the time- dependent 
Hamiltonian given in Eq. (PQ). We can diagonalize Eq. ([1]) as 



H= - 



-{P\t)P{t)-Q\t)Q{t)\ 



(2) 



where the time-dependent annihilation operators are defined as P = Ub + Va, Q = 
-Vb + Ua, and \U\ 2 = [1 - e{t) / sj e{t) 2 + A 2 }/2, \V\ 2 = [1 + t{t)/^e{t) 2 + A 2 }/2. Note 
that \U\ 2 + \ V\ 2 = 1. We assume the system is initialized in the state \N) p \0)q. 

Now suppose we divide the evolution time from —T/2 to T/2 into n time intervals. 
At tk = —T/2 + k5t, where St = T/n and k = 1, • • • ,n, we perform the measurement 
on each boson in the diagonal basis of the temporal Hamiltonian. Suppose at t = tk 
the system is projected to \N)p k \0)Q k . The probability that the system is projected to 
|0> Qfc+1 at t = t k+1 is given by 

= (U k+1 U k + V k+1 V k ) N . (3) 

vSt = 2eo/n and the dimensionless parameter 5r k+ i = 
5t/^t{t k ) 2 + A 2 < 5e/A which can be arbitrary small when n is large enough (i.e. 
the frequency of measurement is high enough). Then we can approximate Eq. ([3]) as 



Probfc. 
We define 8e 



Prob 



fc+i 



A 2 



8(e(t k ) 2 + A 2 ) 



+ o{Ar 



2 > 



N 



> 



N 



The probability of the state being in the ground state at t 
bounded by 



Prob(|A0pj0) Q J>nProb A 



> 



k=l 



1 - 



2A 2 



Nn 



^ ■ (4) 

- T/2 is then 
(5) 



The RHS goes to 1 when n goes to infinity, therefore 
lim Prob(|JV)pJO> Q J = 1. 



(6) 



This result shows that the quantum Zeno effect is also capable of stabilizing the Landau- 
Zener transition by discrete projective measurements. The success of the quantum Zeno 
effect only depends on the measurement frequency n and the ratio e^/A . Stabilization 
of the ground state can always be achieved by choosing n sufficiently large, as long as 
there is a A finitenon-zero minimum energy gap A. 
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3. Continuous model 

Although such projective measurements can in principle stabilize the Landau-Zener 
transition, practically such a sequence is very difficult to implement. The reason for 
this is that knowledge of all the instantaneous eigenstates along the adiabatic sweep is 
required in order to carry out such a procedure. This is in general more difficult than 
finding the solution state, which is just one of the eigenstates. 

To utilize the quantum Zeno effect in a practical quantum computation setting, 
let us now include a system-reservoir coupling. Suppose that we continuously evolve 
the system Hamiltonian and open the system to the environment. The decoherence 
induced by the system-reservoir coupling is mainly in the diagonal basis of the system 
Hamiltonian, assuming that the system-reservoir coupling is much faster than the 
evolution of the system Hamiltonian. The coupling to the reservoir can be described 
as a continuous measurement in the diagonal basis of the instantaneous Hamiltonian 
implementing the quantum Zeno measurements. This passive measurement is much 
more practical since it is achieved without knowing the actual eigenstates of the 
Hamiltonian. 

By evolving the system Hamiltonian according to the scheme of quantum adiabatic 
computation, the continuous measurement induced by the system-reservoir coupling can 
realize the quantum Zeno effect and pin the system in the ground state. This evolves the 
system to the final ground state, which is the solution of the computational problem. 

The decoherence due to the system-reservoir coupling is usually decomposed into 
two types: a <r x -type longitudinal relaxation and a cr z -type transverse relaxation. The 
total master equation for both types of decoherence at zero temperature is written as 



r. 



f t = -i[H,p] + -fV \J-\ p +-±V [J z \ p, (7) 

where H is given in Eq. fl2]), and T>[c]p = —&cp — p&c + 2cp& is the Linblad 
operator. J~ = P^Q is the operator for the a^-type longitudinal relaxation, and P x 
is the associated relaxation rate assumed to be constant. J z = (P^P — Q^Q)/2 is the 
operator for the <r z -type transverse relaxation, and P z is the associated relaxation rate 
also assumed to be constant. We consider the system to be at zero temperature for 
simplicity. 

First, we solve the master Eq. ([7]) only with the cx^-type relaxation, by following 
the method in Ref. [12]. The average of any operator A on the system is given by 

j f (A) = iy/e(t)* + A*([j*,A]) + 

-^-{J + J-A-2J + AJ- +AJ+J-). (8) 

We then derive the equation for J* on the i-th boson as 

j t (J-) = ~ «(*)«#> + <JT» - y (#JT* - 2# W + JfJtJT) 
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+ J+J-j* - 2J+J!JJ + JfJtJj). (9) 
The first bracket on the RHS comes from (gjf), where re(t) = (VU - UV) = 
y2u ^~y 2v2 = 2 (e($+A i ) characterizes the strength of the diadiabatic excitation. The last 
two brackets on the RHS come from the o^'-type relaxation. By using Jf = — | + J^Jf, 
the second bracket is evaluated as —r x (~+ (Jf )). The third bracket is then transformed 
into . J.) _ 2(J? J/)), where (J, • J,) = (Jf J? + Jf J/ + Jfjf) can be 

easily proved to be invariant during the evolution. With the initial condition that 
\<i/>(t = -T/2)) = \N) p \0)q, we find that (J< ■ J^) = 1/4, for any j ^ i. Combining 
equations for all Jf, we obtain 

j t (r) = - K(t)((j + ) + (J-)) - r,(±jv + (j 2 )) 

- y [ I n{n ~ 1} " 2 E E< J * J />1- ( 10 ) 

Using Jf = -i + J+JJ again, E^(^) = (JV - l)/4 - (JV - 1)( Jf) + N{N - 
l)(Ji~ Jr Jj' Jj)., for any «' ^ j. Finally, we use the first-order Hartree approximation [12], 
i.e., (J+JrJ+Ji) « (J+^)(J+Jr) = (1/2+ (J z )/Nf, to obtain 

A similar method is used to derive the equations for (J^) with a^-type relaxation. 
We first obtain the equations for (Jf) on the i-th boson as 

-^{J^-r^JtJTJf)- (12) 

Then we employ a similar approximation, namely, {./,' ./ J~- J'- ) rs (Jf^ Jf)(J^ Jj) and 
(J+J> jf) « {J?J7){Jf) for any i ^ j. The result is 

-r.ov-ixJ*) f^ + ^V (is) 

Second, we solve the master Eq. (J7|) only for cr 2 -type transverse relaxation. Similar 
to the first case, we use 

|(i) = s VWT^([J„i]) + (|i 

-y(JV 2 i + iJ 2 J 2 -2J 2 AT). (14) 
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We can easily derive the equations for the angular momentum for a single boson and 
extend the results to those on the collective operators as 

j t (J z ) = ~ <t) ((J + > + (J-)) , (15) 

j f (J ± }= ±iVW^{J ± ) + Mt)(J Z )- I j{J ± )- (16) 

Equations ([IT]) . ( [13]) . ([15]) and ([IB]) show that if the relaxation rates and F x 
are large enough, (J ) is kept as 0. Therefore the system is almost completely in the 
ground state at all times and (J z ) is close to —N/2, even if the evolution is faster 
than the adiabatic limit. This shows that the Landau-Zener transition can be stabilized 
by the continuous quantum Zeno effect induced by a strong enough system-reservoir 
coupling. Note that a finite minimum bandgap A is also necessary to make the Landau- 
Zener transition successful with the quantum Zeno effect. Without this condition, n(t) 
diverges and (J ) cannot be kept to under any condition. 
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Figure 1. (Color online) The probability of an individual boson being in state a 
during the evolution with a cr^-type longitudinal relaxation for various N and r x . 



Figures [T] and |2] demonstrate the time-dependent evolutions of the Landau-Zener 
transition with a cx x -type longitudinal relaxation and a <r z -type transverse relaxation. In 
both cases, the initial condition is \ip(t = — T/2)) = |JV)p|0)q. Both results demonstrate 
that by increasing the relaxation strength the quantum Zeno effect can successfully 
stabilize the Landau-Zener transition. However, significant differences between the 
two cases exist. First, in Fig. [TJ the presence of many bosons enhances the effective 
relaxation strength as can be seen by examining the dependence with N in Eq. (ITT]) and 
( fl~3l) . Even if T x is very small, a large N makes the system evolution approach to the 
"quantum Zeno effect limit" . On the other hand, for the case of cx^-type relaxation in 
Eqs. (TT5]) and f[T6]) , the transverse relaxation rates are independent of N. Thus many 
bosons will not enhance the final probability of the system being in the ground state. 

Second, if the effective relaxation strength is not large enough, the u^-type 
relaxation can still drive the system back to the ground state by dissipation, for example, 



Enhanced quantum Zeno effect in the bosonic system 



7 



r = 100 

z 

r = 10000 

z 

£ o = 100 
A = 0.3 
T = 100 



't 

-50 -40 -30 -20 -10 10 20 30 40 50 
Time 

Figure 2. (Color online) The probability of an individual boson being in state a 
during the evolution with a <7 z -type transverse relaxation for various T g . 

by emitting phonons to the reservoir. The relaxation is also enhanced by larger N, by 
comparing the dashed line and the dashed-dotted line in Fig. [TJ Cases where the system 
is significantly excited at the anti-crossing point (t = 0) but relaxes to the ground state 
after the evolution finishes may be called the "bosonic final state stimulation limit" . The 
case of Fig. [2] does not have such effect since it only causes dephasing to the system. 
Summarizing the two points, the cr^-type relaxation is advantageous over the <r z -type 
relaxation. With many bosons, both cases can improve the measurement signal-to-noise 
ratio, but the a^-type relaxation will further boost the measurement sensitivity. 

4. Success probability 

We now turn to calculate the success probability of the Landau-Zener transition with 
continuous quantum Zeno measurements. As we have shown that the cr x -type relaxation 
is advantageous, we will focus on this type of relaxation to demonstrate the speedup 
with many bosons. 

In the Landau-Zener transition, the answer state is the up-spin state a. In tge 
many-boson system, the algorithm outputs state a if more than half of the bosons are 
measured to be in state a. Given that p is the probability of an individual boson being 
in state a and the measurement outcome of each boson is either or 1, (state b or a 
respectively), for large N, the ensemble average of the measurement outcome of each 
boson can be approximated by a Gaussian distribution with mean value \x = Np and 
variance a 2 = Np(l — p). As the cut-off is 0.5, the probability that the algorithm 
fails is P e = [1 - erf(^)]/2, and the success probability is 1 — P e . Therefore, given 
P e < 10~ 12 , a typical tolerable error rate in communications, we can find the requirement 
for p and thus simulate the minimum evolution time to achieve this requirement. Signal- 
to-noise ratio can be improved by increasing N, in the sense that the lower bound for 
p is decreased. If N is large enough, the improvement of signal-to-noise ratio is less 
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Figure 3. (Color online) Simulation of minimum time T m i„ to reach P e < 10 12 with 
various N. 

significant. Nevertheless, cr z -type relaxation can further boost p (i.e. the measurement 
sensitivity) by adding more bosons. 

Fig. E] illustrates the speedup of the Landau-Zener transition with increasing 
number of bosons in the system. At large N, the relation approaches a straight line on 
the log- log plot. The fitted result gives a polynomial speedup of around N~ - 81 under 
the condition that u^-type relaxation plays a significant role at zero temperature. Note 
that in order to achieve P e < 1CT 12 , the probability of a boson being in state a does not 
need to be close to 1. Therefore, the algorithm always works in the bosonic final state 
stimulation region where the system relaxes to the ground state via actual dissipation, 
for example, by emitting phonons. Although it is only a polynomial speedup, the 
speedup can be in practice quite large for a large number of bosons. For example, in 
proposals where Bose-Einstein condensates are used as the effective qubits, there can 
be more than 10 4 bosons. 

5. Conclusion 

We have demonstrated that the Landau-Zener transition can be stabilized by both the 
discrete and the continuous quantum Zeno effect. In the continuous model, two types 
of relaxation, the <7 x '-type longitudinal relaxation and the <7 z -type transverse relaxation 
were analyzed. The reason for the stabilization of the ground state is that the coupling to 
a reservoir can give information relating to whether the system has been excited above 
the ground state. For example, in the case of strong cx^-type relaxation, a negative 
measurement result, i.e. a lack of excitations in the reservoir (e.g. phonons), is a 
signal that the system is not being excited. On the other hand, for <r z -type relaxation, 
the measurement readout directly gives information of what state the system is in, 
since it is an effective measurement in the diagonal basis of the temporal Hamiltonian. 
In contrast to previous works on adiabatic quantum computation in open dissipative 
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systems [T3J [EH [E2 US] that have been concerned with the robustness and enhancement 
of the evolution, our work emphasizes the relaxation which can play a more significant 
role than the coherent evolution. 

By introducing the many-boson system quantum Zeno effect with cr^-type 
relaxation, we have found that this can significantly enhance both the signal-to-noise 
ratio and measurement sensitivity. The enhancement is found to be a polynomial 
speedup of the number of bosons in the system. By use of a large number of bosons, 
this introduces the possibility of further accelerating quantum computation beyond the 
standard qubit case. 
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